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In order to understand whether internal stresses in glasses are correlated or randomly distributed, 
we study the organization of atomic bond tensions (normal forces between pairs of atoms). Mea- 
surements of the invariants of the atomic bond tension tensor in simulated 2D and 3D binary 
Lennard-Jones glasses, reveal new and unexpected correlations and provide support for Alexander's 
conjecture about the non-random character of internal stresses in amorphous solids. 

PACS numbers: 61.43.-j 61.43.Fs 



I. INTRODUCTION 

Although the existence of internal stresses in glasses is 
well-established jj] , little is known about their statistical 
properties and about their role in the physics of glasses 
and amorphous solids. Intuitively, one expects that since 
local stresses and strains in a glass are a by-product of 
the structural disorder produced by rapid freezing of a 
liquid (or of incompatible local packing and space-filing 
requirements |2| their distribution should be ran- 
dom. A different view is advocated in ref. , where it is 
argued that the local stresses are not random, and that 
their organization is closely related to the low tempera- 
ture vibrational spectra of glasses. Even though the ques- 
tion whether internal stresses in glasses are randomly or 
are subject to some complex ordering is of fundamental 
importance, to the best of our knowledge this issue was 
not directly addressed by experiments to date. A possi- 
ble reason is that the effects of internal stresses are dis- 
cernible only on atomic length scales, whereas on larger, 
but still microscopic length scales, they average to zero 
and the usual description of a glass as an isotropic solid 
with zero average internal stress, applies. 

Unlike the stress singularities produced by lattice de- 
fects (vacancies, dislocations) in crystalline solids, inter- 
nal stresses in glasses cannot be attributed to strongly 
strained core regions embedded in a weakly strained bulk. 
As a consequence, classical elastic theory methods for the 
calculation of strains and of the resultant stresses |(|, are 
not applicable to glassy systems. The problem has been 
solved for periodic solids subjected to arbitrary homo- 
geneous stress, for which one can extend the continuum 
approach down to atomic scales and calculate the local 
stress acting on a primitive cellQ. Lately, it was shown 
that if the internal stresses in an amorphous solid vary 
sufficiently slowly on microscopic length scales, one can 
obtain a similar (to the crystalline case) expression for 
the local internal stress, with the volume of the prim- 
itive cell replaced by the Voronoi volume Q]. However, 
recent numerical studies of atomic stresses in two and 
three dimensional glassy systems at mechanical equilib- 
rium, have demonstrated that the atomic environments 



vary rapidly even between neighboring atoms at zero 
temperature || |[ [l(J. The breakdown of the contin- 
uum approach implies that the macroscopic concept of 
stress can not be extrapolated down to atomic scales in 
glassy solids. 

In order to bypass these difficulties, in this paper we in- 
troduce the atomic bond tension tensor that accounts for 
the magnitudes and the directions of the normal forces 
(bond tensions in the terminology of ref. Q) acting be- 
tween pairs of atoms with radial pair interaction. Since 
this tensor can not be identified with the stress tensor, 
the original question about the random vs. non-random 
character of the internal stresses, is replaced by that 
concerning the organization of atomic bond tensions in 
glasses. 

In section II we define the atomic bond tension tensor 
and discuss its physical meaning. In section III we de- 
scribe the measurement of this parameter in simulated 
two and three dimensional binary Lennard-Jones Inher- 
ent States of liquids (using the terminology of ref. [pi]) 
and of glasses in mechanical equilibrium. In section IV we 
show that the local organization of bond tensions belong- 
ing to an atom is not random and that the observed cor- 
relations between the eigenvalues of the proposed tensor 
have a physical meaning. In section V, directional corre- 
lations and clustering are addressed, and some parallels 
with granular matter are drawn. These results, which 
demonstrate that the atomic bond tension tensor probes 
simulated glasses in a new way, thus extracting informa- 
tion unavailable by the existing methods, are discussed 
in section VI. 



II. THE ATOMIC BOND TENSION TENSOR 

Consider a classical system of atoms with coordinates 
{Ri}, interacting through the pair potential U({Rij}), 



where Ri 



\Ri — RA is the distance between atoms i 



and j. The tension in the "bond" [ij] between the i—th 
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and the j— th atom is defined as 

dRij 

Notice that in a perfect monatomic crystal, all bond 
tensions vanish in mechanical equilibrium, whereas in 
a multi-component crystal there are Jy ^ 0. Non- 
vanishing Jjj are indications of "mismatches" in the 
bonds. 

Whereas the definition of bond tensions (or, equiva- 
lently, of the normal forces) is strictly local, their dis- 
tribution, P(Jij), depends on the spatial order. P(Jij) 
is related to the radial distribution function g(r) by 
P(Jij)dJij = g(r)dr, where P(Jij)dJij is the probabil- 
ity to find bond tension with values between Jy and 
Jij + dJij, and g(r)dr is the probability of finding a par- 
ticle between r and r + dr, given a particle at the origin. 
Therefore, the variation in nearest neighbor separations, 
manifested in the width of first peak of g(r), leads to a 
broad distribution of bond tensions (the distribution of 
the normal forces was calculated for granular materials 
subjected to external forces ]T^|, Q, and for liquid and 
glassy systems pj). The above connection between the 
distributions P[Jij) and g(r) has several important im- 
plications. If the separations Rij between an atom and 
its neighbors are randomly chosen from the distribution 
<?(r), one does not expect to find correlations between 
the corresponding Jy. If separations between nearest- 
neighbors are nearly equal, the corresponding bond ten- 
sions must be nearly equal as well (the atomic environ- 
ments are either isotropically stretched or isotropically 
compressed). 

Although there are numerous possible realizations of 
the set of interatomic distances {Rij}, their number is 
limited by several constraints. In mechanical equilibrium, 
the atomic positions are constrained by the requirement 
that the net force on every atom in the system vanishes, 
i.e., 

A Y^-Ljl'o (2) 
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where Rij is a unit vector in the direction of Rij. More- 
over, mechanical stability requires that every atom in 
the system is located at its local energy minimum. That 
means that not only the Einstein frequencies of every 
atom in the system but also all the eigenvalues of the 
complete dynamic matrix are positive g . These require- 
ments must be satisfied by Inherent States of a liquid 
(in the terminology of ref. as well as by glasses at 

zero temperature. In the case of granular matter ther- 
mal fluctuations play no role and only the requirement of 
force balance applies at rest Jl5| ], |lq ]. This is an impor- 
tant difference, especially in view of a suggestion due to 
Alexander 0), that the structural disorder in glasses at 
zero temperature should be constrained by the stability 
requirements. 



If the distribution of bond tensions belonging to an 
atom i is not isotropic, the directions associated with the 
different bonds should be taken into consideration. This 
is done by introducing the tensor Ty (with i ^ j) with 
components 

— R i; j R^y (3) 

In order to describe the bond tensions acting on an 
atomic environment, the site of atom i is taken as a ref- 
erence point and summation over all Ty (with j ^ i) is 

performed. The resulting atomic bond tension tensor Xj 
is given by 

X 7 '-'- ( 4 ) 

< — > 

The physical meaning of Tj become apparent when it is 
decomposed to its trace 

•A X J H ( 5 ) 
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and to the traceless part Li , with components 

if = E J « ( k ^K-- d ^) (6) 

where d is the spatial dimensionality of the system and 
b a p is the Kronecker delta function. If the distribution of 
atomic bond tensions (ABT) is isotropic, the bond ten- 
sions acting on an atom i are completely described by the 
trace Ji that can be interpreted as the effective pressure 
acting on the atomic environment jL7j, |hJ. The trace- 
less tensor Li , Eq. [| describes the effect of anisotropy 
on the ABT tensor. To illustrate this point, consider a 
monatomic crystal. At mechanical equilibrium with no 
external stresses, all bond tensions vanish. If the sys- 
tem is under hydrostatic compression, the atomic envi- 
ronments vary isotropically so that Ji ^ 0, while all the 
components L^ vanish. Under pure shear, Jj vanishes 
whereas Lf 13 ^ and we conclude that Li describes the 

local effect of shear. Since Li and Ji are related to lo- 
cal shear and local pressure, respectively, for monatomic 
crystals Tj can be identified (upon suitable normalization 
by the surface area of the unit cell) with the local stress. 
For structurally disordered systems, the variation in sys- 
tem's volume would lead, in general, to variations in the 
local Li . If the system is not sheared as a whole, Li 
measures the weighted effect of the intrinsic anisotropy, 
on the distribution of ABT. 

Although the ABT tensor, Eq. || was constructed 
with the aim to characterize the bond tensions acting on 
an atomic environment, its components are not rotation- 
translation invariant and, unlike bond tensions, depend 
on the choice of the coordinate system. In order to 
obtain a basis-independent (scalar) characterization of 
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bond tensions, one can consider the invariants of the ABT 
tensor (e.g., in 2 dimensions, the trace and the determi- 
nant) . A drawback of such an approach is that among the 
d invariants (d is the dimensionality of space), only the 
trace, eq. |^, has the units of bond tension. An alternative 
approach utilizes the fact that the eigenvalues of the real 
symmetric matrix are invariant under arbitrary trans- 
lation and rotation and, therefore, for every atom i in the 
system, there is a set of d eigenvalues {i,} = (t\, ■ ■ ■ , tf) 
and d eigenvectors {e,} = (e|, • •• , ef) of the ABT ten- 
sor. We stress that although {ti} have the same units as 
the components of the force acting on the atom, there is 
a profound difference between the two; the eigenvalues of 
ABT tensor are uniquely defined scalars that can be com- 
pared for different atoms at different times and external 
conditions (temperature, pressure). They do not vanish 
when the force acting on the atom vanishes, and therefore 
are suitable parameters to characterize atomic environ- 
ments in amorphous systems and in Inherent States of 
liquids. 

The object of the present study is to characterize the 
organization of atomic bond tensions in glassy systems. 
To do so using the ABT tensor, the distribution of the 
sets {ti} should be known for every atom in the system. 
Such a distribution is trivial in monatomic crystals under 
hydrostatic compression, for which all the eigenvalues of 
ABT tensor are degenerate. While in multi-component 
crystalline solids the eigenvalues are, in general, non- 
degenerate, the fact that the atoms are arranged on a 
periodic lattice implies the existence of correlations be- 
tween the eigenvalues. Whether the eigenvalues of the 
ABT tensor in structurally disordered materials are de- 
generate, correlated or randomly distributed, is an open 
question. In the following we address this question for 
two and three dimensional binary Lennard-Jones glasses 
in mechanical equilibrium. 



III. THE SIMULATED SYSTEMS 

The simplest (and perhaps the only) way to calculate 
atomic bond tensions is by numerical simulations. Since 
the ABT tensor is a hitherto unexplored quantity, there 
are many factors to be considered when choosing the pa- 
rameters of the simulated system. Local properties, such 
as the existence of short range structural order, or global 
factors, such as the dimensionality of the system [Q, may 
affect the generality of the results. In order to suppress 
short-range-order effects, we use binary mixtures of two 
types of atoms interacting through Lennard-Jones poten- 
tials (results for monatomic systems with a pair poten- 
tial that leads to local icosahedral ordering were reported 
elsewhere JTof). Since 2D models of glasses offer some 
advantages (in terms of computational load and visual- 
ization options) over the more realistic 3D systems, we 
consider both two and three dimensional systems with 
similar repulsive-attractive pair potentials. 

For 3D glasses we use binary 6-12 Lennard-Jones 80 : 



20 mixtures of atoms of types A and B (denoted 3DbLJ) 
with interaction parameters [fl9| : e 



A A 



1.0, e AB = 

1.5, e BB = 0.5, a AA = 1.0, a AB = 0.8, a BB 
0.88. Particles of type A are majority particles (i.e., 80%). 
The interaction vanishes for r > 2.5a AA . 



For 2D glasses we use binary 6-12 Lennard-Jones 50 : 
50 mixtures (denoted 2DbLJ) where ]T^| 
a AB = 0.5(cr A + o B ) for A, B = 1, 2 with r < 4.5a AB . 
All measurements are given in reduced LJ units (for Ar- 
gon these units correspond to e = 120/cb and a = 3.4 
A). The bond tensions are given in units of e AA /a AA for 
the 3DbLJ systems, and £/a 1 for the 2DbLJ systems. All 
atoms are assumed to have unit mass. 



Liquids and glasses were prepared using standard 
molecular dynamics simulations preformed at a constant 
volume p0[. Periodic boundary conditions were assumed. 
For each system type we considered two densities: p = 
1.2 and p — 1.4 for the 3DbLJ systems, and p = 0.596 
and p — 0.735 for the 2DbLJ systems. Such difference in 
densities corresponds to 5% ( 10%) difference in the aver- 
age nearest neighbor separation for the 3DbLJ (2DbLJ) 
system. The number of atoms N for each 3DbLJ sys- 
tem is 16384. For 3DbLJ system with p = 1.2, sam- 
ples with 4000 atoms were also considered. Comparison 
between ABT tensors calculated for the small and the 
large systems shows that similar results are obtained in 
both cases. For 2DbLJ systems, we used samples with 
N = 4096 atoms. 



Since the binary Lennard-Jones systems used in this 
study do not tend to crystallize jl9| , Q , samples at zero 
temperature can be prepared either by steepest descent 
quench from the liquid state (Inherent States of the liq- 
uid), or by obtaining a well-equilibrated glass at a suf- 
ficiently low temperature (at which all atoms are effec- 
tively trapped in their cages) , and then using steepest de- 
scent quench to yield glassy systems in mechanical equi- 
librium. For dense 2DbLJ systems, an additional route 
of preparation was used: taking samples of low-density 
2DbLJ glasses (at zero temperature), we decreased the 
length of the simulation box and rescaled all the atomic 
positions by 10%. Performing a steepest descent quench 
on the rescaled samples (in order to relax the inter- 
atomic forces) yields dense 2DbLJ systems at mechan- 
ical equilibrium. Although the above preparation routes 
are not equivalent, we did not find statistical differences 
(through comparison of means, ranges and variances) be- 
tween samples prepared by different methods. Also, di- 
rectional correlations related to ABT tensor are similar, 
regardless of the preparation method. Details about the 
preparation and the distribution functions P({U}) are 
summarized elsewhere [JL0|. 
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IV. LOCAL ORGANIZATION OF 
EIGENVALUES OF ATOMIC BOND TENSION 
TENSOR 

Since we are not aware of theoretical models that pre- 
dict correlations between bond tensions in structurally 
disordered systems, we first address the question: Is the 
organization of the d eigenvalues of the ABT tensor of 
an atom random? If the atomic environment is such that 
it can be reproduced by a random choice of d numbers 
from the complete distribution function P({tj}), then the 
directional information contained in the ABT tensor is 
redundant. Observation of non-random organization, on 
the other hand, would confirm that the bond tensions 
acting on an atomic environment are correlated. 

First we show that eigenvalues belonging to the same 
atom are not degenerate at zero temperature. This 
is an interesting point because the non-degeneracy of 
the d eigenvalues of ABT tensor implies that the cor- 
responding d eigenvectors form an orthogonal basis, i.e. 

• ef = 8 a p. Since in structurally disordered sys- 
tems the atomic bond tensions corresponding to any two 
(even nearby) atoms i and j, may be very different, one 
does not know a priori how to relate between the cor- 
responding local coordinate systems {eijand {Bj}. If, 
however, the eigenvalues of ABT tensor are not degener- 
ate for every atom in the system, then the local coordi- 
nate system {e{\ for each atom, can be uniquely writ- 
ten in terms of the external global coordinate system 
{Xi, • • • , Xd\ that we choose to project J.^ on. Calcu- 
lation of the fraction / (in %) of pairs of eigenvalues with 
\tf - tf | < 6 where a ^ (3 (a,/3 = 1, . . . , d), yields that 
for S = 0.5, / w 1%, 0.2% and 2 - 3% for 3DbLJ with 
p = 1.2, 3DbLJ with p = 1.4 and 2DbLJ with p = 0.596, 
respectively. For 2DbLJ with p = 0.735, / < 0.1%. 
Taking 8 = 0.1 yields / < 0.1% for all systems except 
2DbLJ with p = 0.596 for which / w 0.1 - 0.15%. These 
calculations indicate that the degeneracy increases when 
the density decreases. Nevertheless, the actual values of 
/ are very small for all the systems considered, proving 
that the probability of degenerate pairs of eigenvalues is 
negligible. 

Next, we compare between the distribution of the 
TV sets {ti} (for each atom i out of the N atoms in 
the sample) and the distribution of the N sets {t r }, 
where the d numbers t^, ■ ■ ■ ,t~ are randomly chosen 
from the measured distribution P({ti}). For each set 
of d numbers, we choose two numbers according to 
tmin = min^ 1 , • ■ • i max = maxft 1 , • • ■ , t d ), and 

(Mmin _ = n ™.(^ • • • , 4), (Mmax = max (*n ' ' ' ,4), 

respectively |2 1| . Contour plots of the non-normalized 
distributions P{{t min , t max }) and P{{(t r ) min , (U) max })) 
for the systems considered (all at zero temperature) are 
shown in figures 1-4. The solid lines correspond to 
tmin — tmax- For the 3DbLJ glasses, each contour plot 
was calculated using one sample with N = 16384 for 
every system. Similar contour plots were obtained for 
other samples of the system (at the same density). In 



order to improve the resolution, in the case of 2DbLJ 
glass we considered sets with N = 16384 and combined 
four samples for each system (we verified that P({t\ 1 £2}) 
are similar in these sets). Comparing the contour plots 
of 3DbLJ p = 1.2 and p = 1.4, and 2DbLJ p = 0.596 
and p = 0.735 systems (figures l.a, 2. a, 3. a and 4. a, re- 
spectively), with the contour plots of the corresponding 
random sets (figures l.b, 2.b, 3.b, 4.b, respectively), we 
find that the distribution of eigenvalues of the ABT ten- 
sor differs from that of the random sets. In particular, 
the non-degeneracy of eigenvalues {t m i ni t max } is very 
pronounced in the simulated 3DbLJ and dense 2DbLJ 
glasses, but not in the corresponding random sets. We 
conclude that the observed non-degeneracy of eigenval- 
ues belonging to the same atom is the consequence of a 
physical constraint, and not just a trivial consequence of 
the small probability to randomly choose d numbers with 
close values. 

Examining the simulated and random histograms, one 
observes that the peak of P({t m i n , t max }) is considerably 
higher than the peak of P ({(t r ) min , (*r) max })- The differ- 
ences between the most probable "simulated" and "ran- 
dom" combinations of {t m i n , t max \, can be estimated by 
comparing the position and the height of "random" and 
"simulated" peaks. Those values for the systems in fig- 
ures 1-4 are summarized in the table below. The fraction 
/ of atoms with P(t m i n ,t max ) > 2/3h s , where h s is the 
height of the simulated peak, is also presented: 





peak coordinates 


peak height 


f(P > 2/3h s ) 




simulated 


random 


simulated 


random 


simulated 


random 


1 


(-9, 5) 


(-7, 4) 


100 


80 


0.27 


0.06 


2 


(-50, -12) 


(-44, -14) 


75 


50 


0.24 


0.003! 


3 


(1,3) 


(1, 1.5) 


200 


190 


0.36 


0.19 


4 


(-44, -32) 


(-42, -40) 


100 


90 


0.20 


0.13 



Comparison between the peak coordinates indicates that 
the most probable {t m i n , tmax} are similar for the sim- 
ulated and the corresponding random distributions. For 
the 2DbLJ and 3DbLJ with the lower densities, the corre- 
sponding t m i n (tmax) have small negative (positive) val- 
ues, indicating that the average atomic bond tensions are 
small. For the dense systems, the corresponding values 
of tmin and tmax are both negative, reflecting the com- 
pression of typical atomic environments in these systems. 
On the other hand, comparison of the peak heights and 
of the fraction of atoms in the vicinity of peaks, demon- 
strates that the probability for the preferred combina- 
tions of {t m i n ,t max } that was measured for simulated 
systems, can not be reconstructed by a random choice of 
eigenvalues. 

There are also marked differences in the shapes of 
the distributions: the contours of P({t m i n ,t max }) re- 
semble a cone that expands towards smaller values of 
t min , whereas those of P({(t r ) min , (t r ) max }) tend to have 
a more circular form. To quantify the difference be- 
tween the respective histograms, we calculated the Pear- 
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son product moment correlations 

PX ' Y ^{X.-XY/n^^-YY/n 

between t min and t max and between (t r ) min and (U) max - 
The obtained values are (with last significant digit 0.05): 



system 




o(t r ■ r 1 

1 \ min ' max J 


3DbLJ (p = 1.2 ) 


0.7 


0.3 


3DbLJ = 1.4 ) 


0.75 


0.3 


2DbLJ (p = 0.596) 


0.75 


0.5 


2DbLJ (p = 0.735) 


0.7 


0.5 



Although the shapes of the distributions (compare figures 
la and 2a for the 3d case, and figures 3a and 4a for the 
2d case) and the correlations between t m i n and t max are 
quite similar for the densities studied, the dispersion in 
the values of t m i n and t max is much larger for the denser 
systems (notice the different scaling of the axes in the 
above figures). 

The observed high correlations show that if t m i n has 
a small (negative) value then it is highly probable that 
tmax is also small. Vice versa, if t m i n is large we expect 
tmax to be large too. The absence of local environments 
with t m i n <C and t max 3> in the simulated systems, 
suggests that the local organization of the eigenvalues of 
ABT tensor is more isotropic than if they were randomly 
distributed. On the other hand, the previously discussed 
observation that the eigenvalues belonging to the same 
atom are always non-degenerate, implies that this ten- 
dency towards local isotropy of atomic bond tension is 
frustrated, and that the local stresses in the glass state 
can not be completely characterized by isotropic pres- 
sure. 



V. DIRECTIONAL CORRELATIONS AND 
CLUSTERING 

Having shown that the eigenvalues of ABT tensor be- 
longing to the same atom are correlated on average, we 
address a related question of whether atoms with similar 
eigenvalues U are clustered in space. Since for most of 
the atoms the tensions in their bonds are quite small, the 
issue of clustering is mostly relevant for atoms with ex- 
treme (high and low) atomic bond tensions, i.e. for the 
tails of the distribution of P{ti). Using different sam- 
ples of 2DbLJ and 3DbLJ glasses at zero temperature, 
we looked (for each sample) for clusters of atoms with 
tmin (t max ) below (above) some given threshold. We 
did not observe clustering of atoms with high t m ax in 
any of the samples. For dense 2DbLJ glass and for all 
3DbLJ glasses studied, atoms with the lowest (negative) 
values of t min were found to cluster together. 

A simulated realization of a dense 2DbLJ glass in me- 
chanical equilibrium is presented in figure 5 (both types 



of atoms are drawn with the same size). The 33% 
of atoms with the lowest (negative) values of t m i n are 
marked with yellow color and the remaining 67% are col- 
ored in blue, light blue, white, pink and red, correspond- 
ing to increasing values of t m in- As can be seen from the 
figure, the yellow atoms form nearly one dimensional per- 
colating chains. Similar pictures were obtained for other 
samples. The measured percolation threshold for atoms 
with lowest t m in is ~ 30%. This value is lower than 
the theoretical percolation threshold for two dimensional 
lattices p2j. Indeed, for randomly chosen 30% of atoms, 
we found considerably smaller clusters (with less atoms) , 
that did not percolate. 

Repeating the above measurements (with different 
thresholds) for snapshots of dense 2DbLJ glasses at 
higher temperatures, we did not observe such "force 
chains" . On the other hand, running these systems 
at temperatures well within the vibrational regime (i.e. 
when no hopping occurred) , and then performing a steep- 
est descent quench to freeze the atomic positions at local 
energy minima (see details in |l(|), we observed similar 
patterns. This observation suggests that the mechani- 
cal equilibrium condition of zero force on each atom is a 
necessary condition for the formation of percolating force 
chains. 

To further investigate this point, a different prepara- 
tion scheme was used to generate dense 2DbLJ systems. 
Using samples of frozen 2DbLJ glass with p = 0.596, all 
atomic positions and the length of the simulation box 
were rescaled by 10%, then the system was frozen again 
to its new local minima. Comparing between the initial 
(with p = 0.596) and final (with p — 0.735) samples, we 
observed that atoms with smallest t m in are not the same 
in the original and the final states. We also observed that 
such compression alone can not account for the formation 
of force chains. Only after relaxation to a state of lo- 
cal mechanical equilibrium, the resultant sample formed 
"force chains" similar to those previously described. 

A percolating cluster of atoms with lowest values of 
tmin in 3DbLJ glass at zero temperature is shown in fig- 
ure 6 (to improve visualization, only a single cluster is 
plotted). If the threshold for t m i n is chosen to include 
~ 5 — 6% of atoms (we found that with such a threshold 
the system is near to or slightly above percolation), the 
percolating clusters form nearly one dimensional chains 
and loops. Although the threshold values of t m i n are 
different for the compressed and non-compressed 3DbLJ 
glasses (t m i n < —77.5 and t m in < —23, respectively), we 
did not observe qualitative differences in the force chains. 

From the observations described in the previous para- 
graphs, we conclude that both the mechanical equilib- 
rium condition of zero net force on every atom and small 
values of £ m j„, are necessary for the formation of force 
chains. To see why this is so, recall that mechanical 
equilibrium implies that the components of the force on 
an atom along mutually orthogonal directions have to 
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vanish. In particular, for a = 1, . . . , d 

N 

ft Y^-i, ,<;,■(■; o. (7) 

3=1 

where f a is the component of the force acting along the 
eigenvector e". The environment of a typical atom i 
(corresponding to the vicinity of the peak of the dis- 
tribution shown, e.g., in figure la) consists of mildly 
compressed and mildly stretched bonds, and is therefore 
nearly isotropic. In this case, the mechanical equilibrium 
condition, equation |^, can be satisfied by balancing the 
negative and positive tensions in the bonds of the atom 
with all its neighbors. Now, consider an atom with 

N 

(*i)min = min^ Jij (fij ■ if) 2 < 0. 

3=1 

that belongs to the tail of the distribution. Inspection of 
figure la reveals that for t min <C 0, the spread of values 
taken by i max (and, in 3d, by the intermediate eigen- 
value) increases with decreasing t m [ n resulting in increas- 
ingly anisotropic atomic environment. In this case, the 
condition of zero force on atom i can only be satisfied by 
two similarly compressed bonds with atoms j and k, ori- 
ented at nearly opposite directions. In order to balance 
the forces on atom j there should be another nearest- 
neighbor with a highly compressed bond in a nearly op- 
posite direction, and so on, leading to the formation of 
force chains. Note that above arguments imply that force 
chains would also be observed if one considers extremal 
normal forces Jij instead of extremal eigenvalues of the 
ABT tensor ti. This result is in accord with that of ref. 

S- 

The above arguments can also explain, why for two 
and three dimensional binary LJ glasses, similar forma- 
tion of force chains of atoms with high values of t max is 
not likely (although in ref. H weaker force chains were 
observed for pairs of highly stretched bonds) . For these 
systems, the eigenvalues of ABT tensor corresponding to 
a highly stretched atomic environment are of similar mag- 
nitude (see figures la, 2a, 3a and 4a) and force balance is 
achieved by nearly isotropic stretching of bonds between 
nearest neighbor atoms. Such a nearly isotropic organi- 
zation implies that even if atoms with high t max were to 
cluster, the clusters would not take the form of nearly 
one-dimensional chains, but rather that of isotropic lo- 
calized inclusions. These results demonstrate that the 
directional organization of highly stretched atomic en- 
vironments is different from that of highly compressed 
ones. Since t a are related to the internal stresses that are 
present in glassy systems even in mechanical equilibrium, 
our numerical results agree with Alexander's conjecture 
about the different organization of positive and negative 
internal stresses in glasses (the atoms with high values of 
tmax correspond to the " soft inclusions" discussed in ref. 

§)■ 



VI. DISCUSSION AND CONCLUSIONS 

In this paper we defined a new parameter, the atomic 
bond tension (ABT) tensor, and used it to study the 
directional organization of atomic bond tensions. We 
found that in sufficiently dense systems at very low tem- 
perature, the eigenvalues {t±, ■ ■ ■ , to} of ABT tensor are 
correlated i.e. are on the average, similar in sign and 
magnitude. Thus, if one observes that the environment 
of an atom is stretched (compressed) along some direc- 
tion, one can predict with high probability that it would 
be compressed (stretched) along other directions as well. 
Based on this correlation, one may conclude that atomic 
environments in glasses in mechanical equilibrium are 
isotropic, at least in a statistical sense. However, since 
the eigenvalues of the ABT tensor of an atom are strictly 
non-degenerate, we conclude that the tendency towards 
perfect isotropy of atomic environments is frustrated. 

Since for systems with radial pair potential there is 
one-to-one correspondence between the separation iiy 
and the tension Jij in the bond [ij], the observed cor- 
relations between the bond tensions imply also corre- 
lations between the corresponding Rij. Thus, various 
combinations of inter-atomic distances, and especially of 
nearest-neighbor separations in glassy systems, are con- 
strained, and their local organization is not random. The 
observed deviations from randomness are associated with 
directional correlations; if only isotropic measures of dis- 
order such as local volume strain JH, deviation from local 
sphericity Q or local pressure |L7j are considered, such 
correlations are not taken into account and a more ran- 
dom description of dense liquids and glasses results. 

The results of numerical simulations reported in this 
paper provide a direct confirmation of Alexander's con- 
jecture that "the reference state of the glass cannot be 
regarded as an arbitrary random state" (see chapter 16 
in ref. 0]). Qualitatively, correlations among ABT are 
expected from the requirement that the forces on each 
atom of a structurally disordered system in mechanical 
equilibrium are balanced, and that each atom is caged in 
a potential well, i.e., that each atom is subject to a set of 
constraints that include zero net force and a stable atomic 
environment. For systems with a radial pair potential, 
the one-to-one correspondence between the separation of 
two atoms and the tension on their "bond" means that 
achieving balance of forces and local stability simultane- 
ously for every atom in the system, reduces the number 
of degrees of freedom and sets strong constraints on the 
observed ABT. The observed strong correlations between 
the eigenvalues of ABT, in two- and three-dimensional bi- 
nary LJ glasses, imply that when the asymmetry between 
the bonds is strong (e.g., due to two types of atoms), 
the directional organization of bond tensions that would 
comply with force balance and local stability is subject 
to stronger constraints. 

Since the Einstein frequencies of atomic vibration 
around an equilibrium position are determined by the 
derivatives of the inter-particle pair potential, the ob- 
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servation of correlations between the corresponding Rij 
implies that the atomic force constants are not random 
as well. Thus, further quantification of correlations in 
atomic force constants is required before one can incor- 
porate them into lattice-based models which aim to re- 
produce the inhomogeneous network of force constants 
measured for simulated glasses (see, e.g., [^4| ). 

The self-organization of normal forces under various 
external stresses (compression, shear) has been exten- 
sively studied for granular materials 125, 26 27], |2q] , and 
recently also for liquids and glasses ]14|, and for model 
random packing p9| . We studied the self organization 
of bond tensions on an atomic level, and showed that 
external compression has strong qualitative effect on the 
local organization of the atomic bond tensions. In com- 
pressed two and three dimensional binary LJ glasses, 
atoms with the smallest eigenvalues of ABT tensor are 
organized in nearly one dimensional chains and loops. 
Since force chains were observed only for systems at 
zero temperature, we concluded that balance of forces 
on each atom is a necessary requirement for such orga- 
nization. Although this is a customary requirement in 
modeling of force chains in granular materials [ jl2| , to 
the best of our knowledge such a requirement, although 
obvious, was not previously suggested in the context of 
glasses. In glasses with repulsive-attractive pair inter- 
actions, force balance can be achieved both by direc- 
tional organization of stretched and compressed bonds 
in a way that their total effect is zero, or by strong, 



nearly collinear compression of pairs of bonds (similar 
stretching of pairs of bonds is not possible, since it is 
not a stable configuration). For externally compressed 
glassy systems, there are not enough stretched bonds to 
balance the compressed ones, and only self-organization 
into one-dimensional structures (force chains) of highly 
compressed bonds can balance the forces for every atom. 

Finally we would like to comment on several open is- 
sues for future research. We do not yet understand the 
physical origin of the observed conflict between the direc- 
tional correlations and the strict non-degeneracy of the 
eigenvalues of ABT, and of the resulting frustrated ten- 
dency towards local isotropy of atomic environments in 
glasses. While we demonstrated the existence of organi- 
zation of the atomic bond tensions in glasses and in Inher- 
ent States of liquids and quantified the extent to which 
they deviate from a random distribution, a more intuitive 
understanding of the difference between such complex or- 
ganized states and purely random ones, is clearly needed. 
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Figures 

1. a. Histogram of P{t min ,t max ) for 3DbLJ system 
with p = 1.2 at T = (N = 16384). b. Histogram 
p (tmin^max) for thc same 3DbLJ system. 

2. a. Histogram of P(t m i n ,t max ) for dense 3DbLJ 
system with p = 1.4 at T = (N = 16384 par- 



ticles), b. Histogram P(t r mm ,t r max ) for the same 
dense 3DbLJ system. 

3. a. Histogram of P(t m i n ,t max ) for 2DbLJ systems 
with p — 0.596 at T = 0. (Total number of particles 
is 16384.) b. Histogram P {t r m i n ,t r max ) for the same 
2DbLJ systems. 

4. a. Histogram of P(t m i n ,t max ) for dense 2DbLJ sys- 
tems with p — 0.735 at T = 0. (Total number of 
particles is 16384.) b. Histogram P(t r mm ,t r max ) for 
the same dense 2DbLJ system. 

5. t m i n for a sample of dense 2DbLJ system at zero 
temperature. Atoms marked in yellow are those 
with t min < — 56 (this subset contains 33% of atoms 
with lowest t m i n ). 

6. 3DbLJ systems with p = 1.4 - loop-like cluster of 
atoms with t m in < —77.5. 
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